The Star Formation History and the spatial distribution of stellar 
populations in the Ursa Minor Dwarf Spheroidal Galaxy 

Ricardo Carrera 

(Si '■ Institute de Astroffsica de Canarias, E38205-La Laguna, Tenerife, Canary Islands, Spain 

O 
O . 

^ . Antonio Aparicio 

Instituto de Astroffsica de Canarias, E38205-La Laguna, Tenerife, Canary Islands, Spain 
Q>^ . Departamento de Astroffsica, Universidad de La Laguna, E38200 - La Laguna, Tenerife, 

Canary Islands, Spain 



David Martfnez-Delgado 



> 
O 

o 

■ Instituto de Astroffsica de Canarias, E38205-La Laguna, Tenerife, Canary Islands, Spain 

O 
(N 

O ■ Javier Alonso-Garcfa 

Qh' Department of Astronomy, University of Michigan, Ann Arbor, MI 48109 

O 
^. 

! Received ; accepted 

> 



- 2 - 
ABSTRACT 

As a part of a project devoted to the study of the Ursa Minor dSph, the star 
formation history of the galaxy is presented in this paper. The analysis uses 
wide field photometry, encompassing about 1° x 1° (the total covered area being 
0.75 deg^), which samples the galaxy out to its tidal radius. Derivation of the 
SFH has been performed using the synthetic partial model technique. 

The resulting SFH shows that Ursa Minor hosts a predominantly old 
stellar population, with virtually all the stars formed earlier than 10 Gyr ago 
and 90% of them formed earlier than 13 Gyr ago. Nevertheless, Ursa Minor 
color-magnitude diagram shows several stars above the main, old turn-off 
forming a blue- plume (BP). If these stars were genuine, main-sequence stars, 
Ursa Minor would have maintained a low star formation rate extending up 
to ~ 2 Gyr ago. However, several indications (relative amount and spatial 
distribution of BP stars and difficulty to retain processed gas) play against this 
possibihty. In such context, the most rehable hypothesis is that BP stars are 
blue-stragglers originating in the old population, Ursa Minor hence remaining 
the only Milky Way dSph satellite to host a pure old stellar population. A 
marginally significant age gradient is detected, in the sense that stars in outer 
regions are slightly younger, in average. 

The distance of Ursa Minor, has been calculated using the magnitude 
of the horizontal-branch and a calibration based on Hipparcos data of main 
sequence sub-dwarfs. We estimated a distance d = 76 ± 4 Kpc, which is shghtly 
larger than previous estimates. From the RGB color, we estimate a mctallicity 
[Fe/H] — —1.9 ±0.2, in agreement with a previous spectroscopic determination. 
No metaUicity gradients have been detected across the galaxy. 
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Subject headings: galaxies: dwarf — galaxies: fundamental parameters — 
galaxies: individual (Ursa Minor) — galaxies: spheroidal — galaxies: stellar 
content — galaxies: structure 
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1. Introduction 

Dwarf spheroidal (dSph) galaxies are objects of high cosmological significance. In 
generic hierarchical clustering scenarios for galaxy formation, such as cold dark matter 
dominated cosmologies (White & Rees 1978; Blumenthal et al. 1984; Dekel & Silk 1986), 
dwarf galaxies should have formed prior to the epoch of giant galaxy formation and would 
be the building blocks of larger galaxies. The dSphs observed today would be surviving 
objects that have not merged within larger galaxies. Therefore, unveiling their underlying 
structure could provide important clues about process of dwarf galaxy formation that are 
now observed at high red-shifts. Local Group galaxies offer the only opportunity of studying 
their evolution in great detail through the observation of their resolved stellar populations, 
which are fossil records of the star formation history (SFH). 

DSph companions of the Milky Way were historically considered as old systems 
inhabited by globular cluster- like, population II stars (Baade 1963). However, later works 
revealed traces of intermediate- age stars, such as carbon stars (Aaronson & Mould 1980; 
Mould et al. 1982; Progel et al. 1982; Azzopardi, Lequeux & Wester lund 1986) and bright 
AGE stars (Elston & Silva 1992; Freedman 1992; Lee, Freedman & Madorc 1993; Davidge 
1994; but see also Martinez- Delgado & Aparicio 1997). More recent studies, based on the 
analysis of color-magnitude diagrams (GMDs) using synthetic GMDs, have shown that 
dSph are objects with complex and varied SFHs (Mighell 1997; Martfnez-Delgado, Gallart, 
& Aparicio 1999; Gallart et al. 1999a; Hernandez, Gilmore & Valls-Gabaud 2000; Aparicio, 
Carrera & Martmez-Delgado 2001) and the idea that no two dSphs could be identified to 
have similar SFHs has become popular. 

As part of a bigger project devoted to the study of the Ursa Minor dSph galaxy (see 
Martfnez-Delgado et al. 2001, 2002), we present in this paper the results on the SFH of this 
galaxy. Our data sample the galaxy out to 40' from its center, which is of the order of its 
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tidal radius (see Irwin & Hatzidimitriou 1995; Kleyna et al. 1998; Martinez- Delgado et al. 
2001). Ursa Minor is the second closest dwarf spheroidal galaxy satellite of the Milky Way 
{d = 69 Kpc) and was discovered in the Palomar survey by Wilson (1955). Previous works 
show an old system similar, in age and abundance, to the ancient metal-poor Galactic 
globular cluster M92 (Olszewski & Aaronson 1985; Mighell & Burke 1999) and owning a 
well populated, predominantly blue horizontal-branch (HB) (see e.g. Cudworth, Olszewski 
& Schommer 1986; Kleyna et al. 1998). Martinez- Delgado & Aparicio (1998) pointed 
out the possibility that a marginal intermediate- age stellar population could exist in the 
galaxy, associated to the blue- plume present in its CMD. Olszewski & Aaronson (1985) and 
Mighell & Burke (1999) derived distance moduh of (m - M)o = 19.0 ± 0.1 and 19.18 ± 0.12 
respectively, from a sliding fit to the M92 ridge hue. As for the metaUicity, it has been 
recently estimated by Shetrone, Cote, & Sargent (2001), from spectroscopic measurements 
of giant stars, to be [Fe/H] = -1.9 ± 0.11. 

Martmez-Delgado et al. (2001) have shown evidences that Ursa Minor is in a process 
of tidal disruption. This motivates the spatially extended study of the SFH that wc present 
here. Before us, Hernandez et al. (2000) (see also Valls-Gabaud, Hernandez, & Gilmore 
2000) have obtained the SFH for the very central part of Ursa Minor using deep, but small 
field (~ 2f5 x 2^5) HST data. They find a stellar population mainly composed of stars older 
than 10 Gyr. 

The paper is organized as follows. In §2, the observations and photometric reduction 
are presented. In §3, the CMD of Ursa Minor is discussed, followed by the estimation of the 
distance and metaUicity. The rest of the paper is mainly devoted to the derivation of the 
SFH and the discussion of possible evolution scenarios for Ursa Minor, to which §4, 5 and 6 
are devoted. Finally, the main results of the paper are summarized in §7. 
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2. Observations and data reduction 

The observations of Ursa Minor were carried out in B, V, R and / Jolinson-Cousins 
filters with the 2.5 m Isaac Newton Telescope (INT) at Roque de los Muchachos Observatory 
in La Palma and, with the 0.8 m, IAC-80 telescope at Teide Observatory in Tenerife. In 
the INT we used the WFC installed at the prime focus, that holds four 2048 x 4096 pixels 
EEV chips. The scale is 0.33"/pix, which provides a total field of about 35' x 35'. Three 
fields were observed in B and R filters encompassing a region of about 1° x 1°, with a 
total covered area of about 0.75 deg^. Table |I| provides a summary of the observations. 
Figure shows the observed fields. The field containing the center of the galaxy (field 
A) was also observed in V and /. These data will be used to estimate the distance and 
metallicity (see §3). Observations with the IAC-80 telescope were done to obtain the 
standard transformations for these bands. A Thomson 1024 x 1024 chip was used with a 
scale of 0.43"/pix, covering a total area of nearly 7.5' x 7.5'. 

EDITOR: PLACE TABLE |I] HERE. 

EDITOR: PLACE FIGURE |l| HERE. 

The images were processed in the usual way. Bias and fiat-field corrections were done 
with IRAF. DAOPHOT and ALLSTAR (Stetson 1994) were used to obtain the photometry 
of the resolved stars. Transformation into the standard photometric system requires several 
observations of standard star fields in such a way that standards are measured in all four 
chips of the WFC. In practice, we observed three fields from Landolt (1992) in the four 
chips and used them to calculate relative photometric transformations from each chip to the 
central one (chip 4). Then, a larger number of standard star fields were measured in chip 4 
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during the observing run to calculate both nightly atmospheric extinctions and the general 
transformation into the standard Johnson-Cousins system. In total, 170 measurements 
of 38 standards were made during the observing run of May 1998. Summarizing, the 
transformations from each chip to chip 4 are: 



(64-61) = 0.329 + 0.02(5 -i?); (7 = 0.005, (1) 

(r4 - n) = 0.837 - 0.307(5 - R); a = 0.010, (2) 

(64 - 62) = 0.443 - 0.006(5 -R); a ^ 0.007, (3) 

(r4 - ra) = 0.833 - 0.268(5 - R)- cr = 0.008, (4) 

(54-63) = 0.367 + 0.021(5 -i?); a = 0.008, (5) 

(r4 - rg) = 0.739 - 0.233(5 - i?) ; cr = 0.008, (6) 



where sub-indices refer to chips, lower-case letters stand for instrumental magnitudes, 
and capital letters for Johnson-Cousins magnitudes. The a values correspond to the fit 
dispersion at the barycenter of point distribution; hence they put lower limits to the 
zero-point errors. 

Transformations from chip 4 instrumental magnitudes measured at the top of the 
atmosphere into Johnson- Cousins magnitudes are given by: 



(5 - 64) = 24.553 - 0.061(5 - R); a ^ 0.009, 
{R - Ti) = 24.393 + 0.257(5 - i?) ; cr = 0.010, 



(7) 

(8) 



where, as before, lower-case letters refer to extinction-free, instrumental magnitudes and 
capital letters to Johnson-Cousins magnitudes. 
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The central field (field A in Fig. |l]) of Ursa Minor was also observed with the INT 
in V and I in June 1999. This run was not photometric, for which the IAC-80 telescope 
was used in May 2001 to obtain the photometrical calibration. Fifteen standards from 
Landolt (1992) and several secondary standards defined in the central region of Ursa Minor 
(covered by chip 4 in the INT observations) were observed to this aim with the IAC-80. 
The resulting transformation equations for the IAC-80 were: 

(l^-^^iAC8o) = 21.110 + 0.055(^-1); a = 0.009 (9) 
(/-^lAcso) = 21.064 + 0.076(V-/); a = 0.006 (10) 

while the final transformation for the INT observations (chip 4) resulted: 

(K-y^) = 24.13 + 0.04(V-/); a = 0.02, (11) 
(7-24) = 24.21 + 0.02(y-/); (x = 0.02, (12) 

As formerly, lower-case letters refer to instrumental magnitudes and capital letters to 
Johnson-Cousins magnitudes. 

Finally, dispersions of the extinctions for each night are about a = 0.01 for each filter. 
Aperture corrections were obtained using a set of ~ 200 isolated, bright stars in the Ursa 
Minor frames, dispersions being of the order of a = 0.01. Putting all the errors together, 
the total zero-point error of our photometry can be estimated to be about a = 0.02 for all 
bands. 

In addition to these errors, ALLSTAR provides, for each star, with the dispersion, a, 
of the PSF fitting. Generally, these residuals do not reproduce the external errors of the 
photometry (see Aparicio, & Gallart 1995 and Gallart, Aparicio, & Vilchez 1996), but they 
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are an indication of the internal accuracy of the photometry. Figure ^ shows these residuals 
as a function of magnitude for each star. For the final photometric list we selected stars 
with a < 0.25. 

EDITOR: PLACE FIGURE | HERE. 

Artificial-star tests have been performed in the usual way (Stetson & Harris 1988; see 
also Aparicio & Gallart 1995) to obtain completeness factors as a function of magnitude. In 
short, 13 000 artificial stars were added to the B and R images of Ursa Minor. The central 
chip (chip 4) of central field was used. These stars were distributed in a grid covering 
the whole image. To avoid overcrowding effects, they were separated about 8'.'5 from each 
other. Colors and magnitudes chosen for artificial stars were obtained from the synthetic 
CMD (computed as described in Aparicio & Gallart 1995) of an old stellar population with 
similar metallicity as Ursa Minor. The resulting completeness factors are shown in Figure 
||. Completeness picks at "V" ~ 23.5. But about 5% to 10% of stars are lost at any brighter 
magnitude. These stars are those happening to lie close to very bright, saturated stars. 

EDITOR: PLACE FIGURE § HERE. 

3. The color— magnitude diagram 

Figure | shows the [{B - R),"V"] CMD of the inner (a < 18') and outer (18' < a < 40') 
regions of Ursa Minor ("V^" is defined as "V" = {B + R)/2). It must be noted that stars 
associated to Ursa Minor exist at galactocentric distances larger than 40'. We restrict our 
plot to stars within this distance to limit foreground and background contamination (see 
Martfnez-Delgado et al. 2001 for more details, including CMDs for outermost regions). 
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Figure |^ shows a map of the resolved stars with the elhpses hmiting the inner and outer 
regions over-plotted. The galaxy center, the ellipses ellipticity {1 — h/a = 0.44) and the 
position angle (53°) have been chosen after Irwin & Hatzidimitriou (1995). 

EDITOR: PLACE FIGURE | HERE. 

EDITOR: PLACE FIGURE | HERE. 

The two CMDs shown in Fig. ^ are qualitatively similar. The most noticeably features 
are the faint, main turn-off at "V" ~ 23, (i? — -R) ~ 0.8; the blue-plume extending above it 
and reaching 'V ~ 21.2; a well populated, globular cluster-like, narrow red giant branch 
(RGB), and a predominantly blue HB, including several stars in the RR-Lyr instability 
strip. All this suggest an old, low metallicity stellar population and a small metallicity 
dispersion. But the blue-plume requires more attention because it could be made of 
blue-stragglers (BS) but it could also be the trace of an intermediate-age population. It will 
be discussed in 

Figure |^ shows the \y — iy\ and \V — /,/] CMDs of stars in the central chip of field 
A. Although shallower and less populated than the [{B — i?)," V'] CMD, they are better 
suited for the estimation of distance and metallicity and will be used for this purpose in the 
following paragraphs. 

EDITOR: PLACE FIGURE | HERE. 

3.1. The distance 

The distance to Ursa Minor can be calculated from the HB magnitude at the color 
of the RR Lyrae variables. A relation between [Fe/H] and My of RR Lyrae must be 
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assumed to this purpose. However, Demarque et al. (2000) have shown that such relation 
is not universal and depends on the HB morphology. To overcome this problem, we have 
compared Ursa Minor HB with that of several globular clusters with similar metallicity and 
HB morphology as the galaxy. The distance to Ursa Minor can then be obtained from the 
difference in magnitude between cluster and galaxy HBs and the cluster distances. 

For this purpose, we firstly need an estimate of Ursa Minor metallicity and of its HB V 
magnitude. The metallicity has been estimated using the color of the RGB and the Carretta 
& Gratton (1997; CG from here on) and Zinn & West (1984; ZW from here on) scales (see 
§|3.2|). The first one, [Fe/H] = —1.75, will be adopted here, since the metallicities of the 



globular clusters used for comparison are given in the CG scale. 

The Ursa Minor HB V magnitude must be a reddening free value, which we will denote 
V^B^g, obtained approximately for the center of the RR-Lyr instability strip, i.e., at about 
{V — I) = 0.4. The HB, "V" magnitudes must be firstly transformed into usual V Johnson 
magnitudes. This has been done from the same Padua stellar evolution models (Bertelli 
et al. 1994) used in The average difference for the interval 0.2 < {B — R) < 0.6 is 
V —" V" = —0.031 ± 0.002. Then, an empirically defined fiducial HB taken from Rosenberg 
et al. (1999) has been fitted to the Ursa Minor HB. It produces V^^^' = 19.84 ± 0.07. 
Furthermore, using the IR dust maps by Schlegel, Finkbeiner & Davis (1998), a reddening 
of E{B — V) = 0.034 is obtained for Ursa Minor. Assuming the relations by Cardelli, 
Clayton, & Mathis (1989), this corresponds to an extinction Ay = 0.104, from which 
V^UMi ^ 74 ^ 0.07 results. 



EDITOR: PLACE TABLE HERE. 



The clusters selected for the HB comparison are listed in Table ^. They have been 
taken from the lAC-Padua Catalogue of Galactic Globular Clusters (CGGC) (Rosenberg et 
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al. (2000a,b). Column 1 gives the cluster name. Column 2 lists the metallicity in the CG 
scale. The reddening-free distance moduli and HB magnitudes are given in columns 3 and 
4. These values are estimated from main-sequence fitting to Hipparcos subdwarfs (Reid 
1999 and references therein). Although clusters with similar metallicity as Ursa Minor have 
been chosen, small differences are still present. To avoid the effects of these differences, 
the globular cluster HB magnitudes have been transformed to those corresponding to the 
galaxy metallicity using V^HB,o,[Fe/H] = ^hb,o + 0.18([Fe/i^]t/A/i - [Fe/H]gc). The resulting 
magnitudes are given in column 5. Column 6 lists the reddening-free magnitude differences 
between Ursa Minor and cluster HBs: AVhb.o = ^Hefo ^ ^Bo[Fe/H]- Finally, column 7 gives 
the distance moduli obtained for Ursa Minor as (m — M)q ^' = (m — M)^^ + Ahb,o- 

The average of values listed in table |^, column 7 is (m — M)o = 19.4 ±0.1, corresponding 
to (i = 76 ± 4 Kpc. This value is larger than the value (m — M)o = 19.24 ± 0.24 given 
by Nemec et al. (1988) using RR Lyrae periods. Comparing Ursa Minor and M92 fiducial 
sequences, Miguell & Burke (1999) and Olzewski & Aaronson (1985) derived distance 
moduli of (m — M)o = 19.18 ± 0.12 and (m — M)o = 19.0 ± 0.1, respectively, also smaller 
than our result. 

3.2. The metallicity 

The metallicity of Ursa Minor has been measured spectroscopically by Shetrone et 
al. (2001), who provided a mean value of [Fe/H] = —1.90 ± 0.11. Previously Olzewski 
& Aaronson (1985) photometrically estimated it to be about [Fe/H] = —2.2, assuming 
that Ursa Minor has the same metallicity and age as M92. Here, we will provide a new 
photometric estimate, based on the color of the RGB. This is a frequently used estimator 
because photometrically observing and measuring the RGB is a relatively simple task. 
However, it is clear that this method can not compete with spectroscopic determinations 
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and recently, Gallart et al. (2001) have shown that these RGB-estimates can be wrong by 
far. In any case, the metalhcity we will derive will, at least, be useful for comparison with 
results for other galaxies existing in the literature obtained with the same method. 

To estimate the metallicity, we have used the relations provided by Saviane et al. 
(2000), based on the position of the RGB of several globular clusters of the lAC-Padua 
CGGC (Rosenberg et al. 2000a,b). We will use the calibration based in the {V — J)/=_3 
parameter; i.e., the {V — I) color index of the RGB at Mj = —3. From figure ^ we 
estimate {V - /)/=_3 = 1.17 ± 0.03. The resulting metallicity is [Fe/H] = -1.75 ± 0.20 
if the CG metallicity scale is used, or [Fe/H] = —2.01 ± 0.14 if the ZW metallicity scale 
is used. Considering both values, we adopt [Fe/H] = —1.9 ± 0.2, which agrees with the 
spectroscopic value. 

The color of the RGB can be also used to test whether galactocentric metallicity 
gradients exist in Ursa Minor, under the reasonable assumption that age gradients, if any, 
have negligible effects on the RGB color. To this purpose, the galaxy has been divided in 
several concentric elliptical annuli and the [{B — R),R] CMDs have been plotted for them. 
The BR plane is used now because of the larger spatial coverage made in the observations 
with these filters. Hyperbolic polynomial fits to points in the interval —1.5 < R < —3.2 were 
used to obtain a fiducial RGB for each CMD. The (5 - R) at Mr = -2.5 {{B - R)r=-2.5) 
has then been used as an estimator of the metallicity (Aparicio et al. 2001). It was 
measured for each fit and the results are plotted in Figure |^. The error bars show the 
internal dispersions of points about the polynomial fits. Inspection of this figure reveals 
that no metallicity gradient can be deduced to exist in Ursa Minor. 



EDITOR: PLACE FIGURE HERE. 
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4. The Star Formation History 

The SFH of a galaxy can be derived in detail from a deep CMD through comparison 
with synthetic CMDs (Aparicio 2001). To obtain the SFH of Ursa Minor, we have used 
here the partial mode/ method, introduced by Aparicio, Gallart, & Bertelli (1997). We have 
closely followed the prescriptions and criteria used by Aparicio et al. (2001), where this 
method was applied to derive the SFH of the Draco dSph. Details on the computation of 
synthetic CMDs and on the SFH computation method itself can be found in Gallart (1999) 
as well as in Aparicio et al. (1997), as well as in Aparicio (2001) where different methods 
are reviewed. 

The SFH can be considered a function of several variables (time, stellar mass, etc). A 
typical approach is assuming it composed by the star formation rate (SFR) function, iplt) 
and the chemical enrichment law (GEL), Z{t). But the initial mass function (IMF), 0(m) 
and a function accounting for the fraction and mass distribution of binary stars, (3{f,q) 
are also relevant. The most powerful methods now in use are able to solve for ip{t) and 
Z{t) simultaneously and treat 0(m) and l3{f,q) as inputs (see Aparicio 2001 and references 
therein). In our case, to simplify and having Ursa Minor low metallicity and low metallicity 
dispersion, we have fixed Z{t) and have solved for ip{t) alone. We have adopted the Kroupa, 
Tout, & Gilmore (1993) IMF and have neglected the contribution of binary stars. 

The metallicity given in j j3.2| , [Fe/H] = —1.9 ±0.11 (Shetrone et al. 2001), corresponds 
to Z = (2.5 ± 0.5) X 10~^. We have adopted a constant metallicity law with this value, 
but allowing some dispersion, to compute the synthetic stars. According to it, synthetic 
stars have metallicities randomly distributed in the interval 2 x 10^^ < Z < 3 x 10~^, 
independently of age. This produces a RGB compatible in color and wideness with Ursa 
Minor's. 

The only remaining function to complete the SFH of Ursa Minor is the SFR, ipit)- 



- 15 - 



It has been solved using the aforementioned partial model method. To apply it, a single 
synthetic CMD is computed, embracing the full possible age interval (e.g. from 15 Gyr to 
Gyr). This is divided in several partial models, each containing a stellar population with 
ages in a narrow interval. The Padua stellar evolution library (Bertelli et al. 1994) is used 
to this purpose. The basic idea is that any SFR can be simulated as a linear combination 
of partial models. 

The solution for the SFR must have associated to it a CMD compatible with the 
observed one. To find such solution, a number of boxes are defined in the CMDs (observed 
and partial models) and the number of stars are counted inside them. The boxes are defined 
in such a way they sample stellar evolutionary phases providing information about different 
ages. Let us call the number of stars in box j in the observed CMD and the number 
of stars in box j of partial model i (the partial model covering the i-th age interval). Using 
this, the distribution of stars in boxes of an arbitrary SFR can be obtained from a linear 
combination of the Nj^, by 

N^^k^a^N^. (13) 

i 

The corresponding SFR can be written as 

i 

where ai are the linear combination coefficients, is a scaling constant, and Ai{t) = 1 if t 
is inside the interval corresponding to partial model i, and Aj(i) = otherwise. Finally, the 
tlj{t) having the best compatibility with the data can be obtained by a fitting of to 
N°, the ai coefficients being the free parameters. 

In practice, several choices of partial model distributions have been tried. However, 
choices involving time resolutions for old stars of the order of 1 Gyr are fluctuating in 
a similar way as shown by Olsen (1999), indicating that data do not contain enough 
information for such small time sampling. On the other hand, the SFR at intermediate 
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ages is so low (if any) that only an average estimate makes sense. In summary, after several 
trials, the adopted solution is based on a three partial model decomposition corresponding 
to the age intervals 2-10 Gyr, 10-13 Gyr and 13-15 Gyr. The lower limit (2 Gyr) is the age 
of an isochrone having the turn-off point at "V" ~ 1.5, corresponding to the upper point of 
the Ursa Minor blue plume. 

Moreover, 11 boxes have been defined in the CMD, as shown in Fig. to characterize 
the distribution of stars. Box 1 samples the old MS population. Boxes 2 to 5 sample the 
BP population. In this analysis we will consider two possibilities in turn, namely that 
all the BP stars are genuine MS stars (hence of intermediate-age) and that all them are 
BS (hence old stars). We will discuss more about this in Box 6 to 8 cover the HB, 
while box 9 samples the region corresponding to intermediate mass, core He-burning stars. 
Here, we will assume that age is the only relevant second parameter to determine the HB 
morphology. However, it must be kept in mind that the nature of the second parameter is 
not fully confirmed and that the blue-ward extension of the HB depends on stellar evolution 
parameters that are not well controlled. This could affect the morphology of the SFR for 
the interval 10-15 Gyr, although the effects on the SFR integral for that period would be 
very limited. 

Box 10 samples the RGB, populated by old and intermediate-age stars born in fact 
in the full, 2-15 Gyr interval considered here. For this reason, it gives no information on 
age resolution, but it provides a strong constraint to normalize the full SFR: independently 
of the age distribution, the integrated SFR for the old and intermediate-age stars must 
be compatible with it. Beside a few AGBs in the same age interval as the RGBs, box 10 
includes a number of foreground stars. This number is estimated from box 11 and simply 
subtracted from the counts in box 10. 



EDITOR: PLACE FIGURE | HERE. 
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The solution for the SFR is found in a global way, considering the number of stars in all 
boxes. In practice, it is not reasonable to look for the best solution (the linear coefficients 
best reproducing the distribution of stars in the boxes of the observational CMD), but 
to search all the solutions providing a stellar distribution compatible with the observed 
one within some interval. The criterion used here is to adopt as good solutions all those 
producing chil (the reduced x^) within x^best + 1- 

EDITOR: PLACE FIGURE | HERE. 

The large spatial coverage of our data offers the opportunity of comparing the stellar 
population in different regions of the galaxy. To do so, we have computed ip{t) in the inner 
(a < 18') and in the outer (18' < a < 40') regions defined in These values approximately 
correspond to the core and tidal radii (e.g. Irwin & Hatzidimitriou 1995 find Vc = 15'8, 
rt = 50'.6 and Kleyna et al. 1998 obtain = 34'). The results are shown in Fig. ^. Error 
bars represent the dispersion of the accepted solutions. 

As in the case of the Draco dSph (Aparicio et al. 2001), the solutions are very similar 
in both regions. However, a marginally significant difference could exist in the fraction of 
stars in the 10 — 13 Gyr interval relative to the 13 — 15 Gyr one, in the sense that the former 
could be relatively more frequent in the outer region of the galaxy. Since the difference is 
at the error bar level we do not further discuss on it. 

On the other hand, the solution is largely dominated by very old stars. Some 90% of 
the stars were formed in an initial, main burst lasting from 15 to 13 Gyr ago and at least 
95% of them were formed before 10 Gyr ago. After that, two possibilities arise. If BP stars 
are normal MS stars, then star formation has gone on at a low rate until recent epochs. 
Alternatively, if BP stars are BS, the star formation has been zero or negligible after 10 
Gyr ago. For the reasons discussed below (§^ we adopt the second possibility as the most 
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reliable, which makes Ursa Minor the only Milky Way satellite lacking an intermediate age 
stellar population. 

The obtained SFR has been used to calculate the total mass in stars and stellar 
remnants (M^) within 18' and within 40'. Both values are quoted in Table ^ and the second 
is used to calculate the dark matter fraction, k. 

The SFR has been also computed by Hernandez et al. (2000) (see also Valls-Gabaud et 
al. 2000) for the central region of Ursa Minor using HST data. They obtained a solution 
qualitatively very similar to ours: a major star formation event at an early epoch forming 
all or almost all the stars in the galaxy. However, Hernandez et al. (2000) find a maximum 
value of 6 X lO^^MQyr^^ at about 14 Gyr and for ages < 10 Gyr while Valls-Gabaud et 
al. (2000), 3 X 10~^MQyr~^, also at ~ 14 Gyr. These authors calibrate their SFR scale by 
normalization to account for the total luminosity of Ursa Minor. To compare with them we 
have to sum the solutions for our inner and outer regions. We obtain a maximum SFR of 
2.7 X lO~^M0yr~^ for the period 15-13 Gyr, quite similar to the Valls-Gabaud et al. (2000) 
value. It is important to stress the compatibility of both solutions, which is more noticeable 
considering the different approaches used by these authors and us. 



5. On the nature of blue stars 

To derive the SFH of Ursa Minor, we have considered in §^ two possibilities about 
the nature of BP stars, namely that they are normal, intermediate-age, MS stars or BS. 
Although it is probably impossible to conclusively ascertain what is the right possibility, a 
few tests can be done to check what of them is the most reliable. We will discuss in turn 
the following points, that should cast light on the problem: (i) the amount of stars in the 
RC region of the CMD; (ii) the relative number of BP stars; (iii) their spatial distribution. 
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and (iv) the availability of gas in Ursa Minor. 

5.1. Stars in the RC region 

If BP stars would be genuine, intermediate-age, MS stars, they should have a 
counterpart in the red-clump region of the CMD, sampled by box 9 (Figure |^). However, 
this criterion is not sensitive enough: only 2-3 stars of this kind are expected to be found in 
this box, while some 20 to 22 old, AGB stars evolving upwards from the HE are expected to 
lie in the same box. This must be added to the fact that some foreground and background 
contamination is expected to affect this box, making the star counts in it more uncertain. 

5.2. The relative number of blue stars 

The number of BS stars relative to a reference population, namely the HE stars, is 
used in the study of the frequency and distribution of BS in globular clusters and field 
stars. We have computed this number for Ursa Minor. Only the inner 18' have been used 
to limit the effects of background and foreground contamination. The resulting number of 
BS relative to HB stars is -F^^ = 1.8. This can be compared with the results obtained for 
globular clusters by Piotto et al. (2002) and for Milky Way halo field, blue, metal-poor stars 
(Preston & Sneden 2000). The former obtain values in the range 0.1 < F#/ < 1.0 for most 
globular clusters in their sample, larger values corresponding to less concentrated clusters. 
The only two exceptions have between 2 and 3. Preston & Sneden (2000) obtain 
F^^ = 4.4. This scenario points to the BS possibly originating in close binary stars, that 
would be more efficiently destroyed in high density globular clusters. This would account 
for the higher values of F^^ in lower density environments. Indeed, the value found here 
for Ursa Minor, which is intermediate between those of globular clusters and that of the 
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Milky Way halo field, would be compatible with this scenario and would agree with the BP 
stars in Ursa Minor being BS stars. 



5.3. The spatial distribution of blue stars 

A further test of the nature of the BP stars is based on their radial distribution. The 
gas coming from old, dying stars would have a radial distribution similar to that of these 
stars. But the fact that a critical density of gas is required to form new stars implies that 
these should preferentially form in the inner part of the galaxy (see, for example, the cases 



of Phoenix and NGC 185 in Martinez-Delgado et al. 1999a, b, respectively). Figure |T0 
shows that the distribution of as a function of galactocentric distance in Ursa Minor is 
fiat, indicating that the BP population is strongly related to the old stellar population and 
that it is very likely formed by BS stars. 



EDITOR: PLACE FIGURE 10 HERE. 



5.4. Gas and intermediate-age stars in Ursa Minor 

Producing an intermediate-age population implies a mechanism to generate and 
conserve gas. Deep observations have failed to detect gas in Ursa Minor (Young 1999, 2000), 
even for large distances from the center of the galaxy (Blitz & Robishaw 2000). However, 
the gas could come, at a low rate, from early generation, dying stars. This could eventually 
allow small, short and recursive star formation bursts that, smoothed and averaged, would 
account for a low rate star formation at intermediate ages. The rate at which this material 
is returned to the ISM can be estimated from the integration of the SFR through 
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rrtiu 

^(t) — / [m — p{'m)](j){m)'ilj[t — T{m)]dm (15) 

J mi 

where mi and m„ are lower and upper integration limits for the stellar mass, p{m) is the 
mass of the stellar remnant of a star of initial mass m, and r(m) is the life time of a star 
of mass m. Note that we are using time increasing toward the past, with the present-day 
value being 0. We assume that all the gas returned to the interstellar medium from SNe 
explosions escapes the galaxy and we will estimate only the gas coming from stellar winds 
of intermediate and low mass stars in the inner Ursa Minor region, where density is larger. 
From a linear fit to Vassiliadis & Wood's (1993) results for stars of metallicity Z = 0.004 
(the lowest they compute) in the mass interval from 0.89 Mq to 5.0 M©, we obtain 
p(m) = 0.489 + 0.093 x rrii, where rrii is the initial stellar mass. Using the SFR obtained 
in for the inner region (a < 18') and for the case in which BP stars are assumed to be 
intermediate-age MS stars; the Kroupa et al. 's (1993) IMF; mi = O.89M0, and m„ = 5.0, 
we obtain in average R{t) = 1.3 x 10^^ MQyr"^ for 2 Gyr< t < 13 Gyr (inner 18'). This 
value is five times the averaged SFR obtained in §|| for that period, V^i3_2 = 2.5 x 10^^ 
Meyr-i. 

However, in a dwarf galaxy such as Ursa Minor, gas is expected to be repeatedly 
removed through SNe explosions. Van den Bergh & Tammann (1991) estimated that, 
for a galaxy of the luminosity of Ursa Minor, the SNIa rate is one per 10^ years. In 
this interval Ursa Minor (inner 18') would accumulate some 130 Mq in all of gas from 
intermediate (if any) and low mass dying stars. Under these conditions, it seems really 
difficult to accumulate gas enough in any limited region of the galaxy to allow stars to form 
at intermediate ages. 

To see this more clearly it should be enough considering that, if gas from intermediate 
and low mass dying stars would distribute across the galaxy as the stars do, the central gas 
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density would reach a value of some 4.4 x 10^^ MqPc"^ after 10 Gyr if not removed by SNe 
explosions. This value corresponds to a particle density of 5 x 10^^ cm^^, much smaller 
than values about 1000 cm~'^, which are typical lower thresholds for star forming regions 
(Shu, Adams, & Lizano 1987). 

6. The SFH scenario in Ursa Minor 

In summary, the tests discussed in §^ support the idea that BP stars in the Ursa Minor 
CMD are BS and that not a significant intermediate-age stellar population is present in the 
galaxy. Consequently, the overall SFH for Ursa Minor can be sketched in the following way. 
In a first phase, lasting from about 15 to 13 Gyr ago with a further extension up to 10 Gyr 
ago, most or all the stars now populating the galaxy were formed. Unused gas would likely 
have been ejected or swept in the subsequent SNe explosions. In a second phase, starting 
upon the end of the first one, dying stars injected gas into the interstellar medium. Some 
of it could have been used to form new stars but, most likely, SNIa explosions would have 
avoided the accumulation of enough gas at any moment of this phase, preventing any star 
formation at intermediate ages. In such a case, blue-plume stars in Ursa Minor would be 
BS. The best solution for the SFH is shown by the full line in Figure ^ Alternatively, the 
maximum SFR compatible with the BP stars being genuine, intermediate- age, MS stars is 
shown by a dotted line in the same figure. 

It must be mentioned that, in our analysis of Draco (Aparicio et al. 2001), we assumed 
that the BP was formed by genuine MS stars. The 2-3 Gyr old burst found in the SFH of 
Draco depends on this assumption and, for this reason, as quoted in Aparicio et al. (2001), 
the value found for of the SFR in that burst is to be considered an upper limit. But the 
burst also relies on the RC and subgiant populations present in Draco's CMD, for which 
the existence of an intermediate-age stellar population in Draco should be considered real. 
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Summarizing, although apparently similar at first glance, Ursa Minor and Draco differ in 
two important properties: their dynamical structures, with Ursa Minor showing a tidal 
tail (Martnez-Delgado et al. 2001) but Draco lacking it (Aparicio et al. 2001) and Draco 
showing an intermediate- age stellar population which is not present in Ursa Minor. 

7. Conclusions 

As a part of our project devoted to the study of the Ursa Minor dSph, the SFH of the 
galaxy is presented in this paper. The analysis uses wide field photometry, encompassing 
about 1° X 1° (total covered area about 0.75 deg^) and the synthetic CMD technique for the 
derivation of the SFH. 

The resulting SFH shows that Ursa Minor hosts a predominantly old stellar population. 
Very likely, virtually all the stars were formed before 10 Gyr ago, and 90% of them 
were formed before 13 Gyr ago. This picture shows Ursa Minor as the Milky Way dSph 
satellite most resembling the original hypothesis of Baade, namely that dSphs were pure 
old, globular cluster-like systems. Nevertheless, Ursa Minor CMD shows a well populated 
blue-plume above the main, old, MS turn-off. If these stars were genuine, MS stars, Ursa 
Minor would have maintained a low star formation rate extending up to ~ 2 Gyr ago. 
However, (i) the relative number of BP to HB stars; (ii) the spatial distribution of BP stars, 
and (iii) the gas availability play against such possibility and strongly favors the BP stars 
being a BS population. In this context, Ursa Minor would remain the only dSph, Milky 
Way satellite hosting a pure old stellar population. 

The SFH study has been done for two regions of the galaxy, namely the one within 
the core radius and the one between this and the tidal radius. No significant differences 
have been found between the stellar populations in both regions except for the fact that the 
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fraction of stars found in the 10 — 13 Gyr interval is marginally larger. 

Besides the SFH analysis, the distance of Ursa Minor has been also derived from 
comparison of the galaxy HB magnitude with those of some globular clusters. These were 
chosen to have the same metallicity and HB morphology than the galaxy. We obtained a 
distance of c? = 76 ± 4 Kpc. The metallicity was also estimated from the RGB color to be 
[Fe/H] — —1.9 ± 0.25, in agreement with the spectroscopic value by Shetrone et al. (2001). 
No metaUicity gradient is detected along the galaxy. 

We are indebted to the anonymous referee of this paper for several comments and, in 
particular, for suggesting a further discussion of the blue-stars problem. That point has 
been largely improved thanks to enhghtening comments from Drs. C. Gallart, G. Piotto, 
G. Preston and J. M. Torrelles. 

This article is based on observations made with the 2.5 m Isaac Newton Telescope 
operated on the island of La Palma by the ING in the Spanish Observatorio del Roque de Los 
Muchachos, and with the lAC 80 telescope operated by lAC in the Observatorio del Teide 
of the Instituto de Astrofisica de Canarias. This research has made use of the NASA/IPAC 
Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology, under contract with the National Aeronautics and Space 
Administration. Javier Alonso-Gracia gratefully thanks to the lAC for their hospitality in 
La Laguna. This research has made use of the Digitized Sky Survey, produced at the Space 
Telescope Science Institute at Baltimore under U.S. grant NAGW-2166. This research has 
been supported by the Instituto de Astrofisica de Canarias (grant P3/94), the DGESIC 
of the Kingdom of Spain (grant PI97-1438-C02-01), and the DGUI of the autonomous 
government of the Canary Islands (grant PI1999/008). 
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Table 1. Journal of observations 



Date Ursa Minor Field Time (UT) Filter Exp. time (s) 



98.05.27 


A 


22:39 


B 


900 


98.05.27 


A 


22:57 


B 


900 


98.05.27 


A 

A 


23:15 


R 


900 


98.05.27 


A 


23:33 


R 


900 


98.05.27 


B 


23:52 


B 


900 


98.05.27 


B 


00:11 


B 


900 


98.05.27 


B 


00:29 


R 


900 


98.05.27 


B 


00:48 


R 


900 


98.05.28 


A 


21:28 


B 


30 


98.05.28 


A 

A 


21:32 


R 


20 


98.05.28 


B 


21:51 


R 


20 


98.05.28 


B 


21:54 


B 


30 


98.05.28 


C 


22:10 


B 


30 


98.05.28 


C 


22:15 


R 


20 


98.05.28 


c 


22:46 


B 


900 


98.05.28 


c 


23:04 


B 


900 


98.05.28 


c 


23:22 


R 


900 


98.05.28 


c 


23:40 


R 


900 


99.06.15 


A 


22:59 


V 


10 


99.06.15 


A 


23:02 


V 


10 


99.06.15 


A 


23:06 


V 


600 


99.06.15 


A 


23:20 


I 


10 


99.06.15 


A 


23:23 


I 


600 
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Table 1 — Continued 



* taken at the IAC-80 telescope at Teide Observatory. 


Date 


Ursa Minor Field 


Time (UT) 


Filter 


Exp. time (s) 


99.06.15 


A 


23:37 


B 


120 


99.06.15 


A 


23:43 


R 


120 


01.05.13* 


A 


03:23 


V 


900 


01.05.13* 


A 


03:30 


I 


900 



Table 2. Clusters and parameters used to estimate the distance. 



Cluster Metalhcity (m - M) 



NGC 4590 
NGC 6341 
NGC 6397 
NGC 7078 



-2.00 ±0.03 
-2.16 ±0.02 
-1.76 ±0.03 
-2.02 ±0.04 



15.23 ±0.16 
14.82 ±0.16 

12.24 ±0.1 
15.30 ±0.14 



T/gC 



15.63 ±0.1 
15.14 ±0.1 
12.39 ±0.1 
15.62 ±0.05 



•^HBJFe/H] 

15.67 ±0.13 
15.21 ±0.13 
12.39 ±0.13 
15.67 ±0.1 



AVhb,o 



4.10 ±0.15 
4.56 ±0.15 
7.38 ±0.15 
4.10 ±0.12 



{m-M) 



19.33 ±0.2 
19.38 ±0.2 
19.62 ±0.2 
19.40 ±0.2 



-31 - 



Table 3. Morphological and integrated parameters of Ursa Minor 



Parameter 


Value 


Reference 




15*^09^.2 


(3) 




67°12:9 


(3) 


IJ,"V"fl 


25.5 ±0.3 


(3) 




15f8± lf2 


(3) 


d 


70 ± 4 kpc 


(1) 


[Fe/H] 


-1.9 ±0.11 


(4) (1) 


Vt 


10.3 ±0.4 


(2) 


Mv 


-8.9 


(1) 


Lv (Le) 


3 X 10^ 


(1) 


MvT {Mq) 


23 X 10^ 


(1) 


Mvt/Lv (Mq/Lq) 


77 


(1) 


M,(a < 18') (Me) 


2.25 X 10^ 


(1) 


M.(a < 40') (M©) 


3.6 X 10^ 


(1) 


K = 1 - (M^r < 40')/Mvt) 


0.98 


(1) 



References. — (1) This paper; (2) Caldwell et al. (1992); (3) Irwin & 
Hatzidimitriou (1995); (4) Shetrone et al. (2001) 
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Fig. 1. — Digitized Sky Survey image of the Ursa Minor region. The observed fields are 
over-plotted. The galaxy is centered about the center of Field A. North is up, East is left. 




Fig. 2. — Residuals of the PSF-to-star fitting as a function of magnitude provided by 
ALLSTAR for the stars resolved in the Ursa Minor frames. 
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Fig. 3. — The completeness factor as a function of "V" for the central chip of field A (see 
Fig. I). "1/" is defined as "1/" = {B + R)/2. 
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Fig. 4. — The CMD of the inner and outer regions Ursa Minor. Upper paneh stars within 
a < 18'. Lower panel: Stars within 18' < a < 40'. "V" is defined as "V" = (S + R)/2. 
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Fig. 5. — Spatial distribution of the resolved stars together with the eUipses used to spatially 
divide the galaxy. EUipses semi-major axis are a — 18' and a — 40'. 
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Fig. 6. — [{V — I),V] and [{V — /),/] CMDs corresponding to the central chip of the WFC 
in Field A (see Fig. 0). 
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Fig. 7. — {B — -R)r=-2.5, the color of the RGB measured at Mr — —2.5, of Ursa Minor stars 
for eUiptical annuh of increasing semimajor axis a. 
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Fig. 9. — The SFR as a function of time for the inner (a < 18') and outer (18' < a < 40') 
regions of Ursa Minor. Thick, sohd hues show the solutions obtained assuming that the BP 
stars are BS. Thin, dotted hues, show the solutions obtained assuming that the BP stars are 
genuine MS stars. Error bars represent the solution dispersions for each age interval. 
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Fig. 10. — The radial distribution of the BP stars relative to HB. Boxes 2 to 5 and 6 to 8 
(as defined in Figure H) have been respectively used to sample the BP and HB stars. 



